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ABSTRACT 


The  collisionless  drift  wave ,  with  a  component  of  the 
propagation  vector  parallel  to  the  magnetic  field,  is  examined 
in  a  nonunif ormly  rotating  cylindrical  plasma.  Terms  involving 
the  ratio  of  the  Alfven  velocity  (C  )  to  the  velocity  of  light  (C) 

A 

are  retained  in  the  analysis  allowing  plasmas  in  different 
density  regimes  to  be  investigated.  The  differential  equation 
describing  the  drift  mode  is  solved  numerically.  Approximate 
analytic  solutions  are  also  given  using  the  WKB  method  and  a 
Weber-Hermite  approximation. 

It  is  found  that  the  velocity  shear  resulting  from  non- 
uniform  electric  and  diamagnetic  rotations  can  provide  the 
dominant  contribution  to  the  growth  rate  of  the  wave.  The  de¬ 
stabilizing  effect  of  nonuniform  rotations  is  further  enhanced  in 

2 

lower  density  plasmas  where  terms  involving  (C^/C)  are  no  longer 
negligible. 
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CHAPTER  I 
INTRODUCTION 

Many  phenomena  which  occur  in  a  confined  laboratory  plasma  cannot 
be  adequately  described  by  using  the  infinite  homogeneous  model.  The 
effects  of  temperature  gradients,  density  gradients,  finite  lengths, 
boundary  conditions  and  external  field  configurations  become  significant 
in  the  laboratory  and  thus  can  produce  a  whole  new  host  of  waves  and 
instabilities.  It  is  worthwhile  digressing  for  a  moment  as  to  the  exact 
meaning  of  an  instability.  Since  the  Hamiltonian  of  a  plasma  does  not 
contain  magnetic  field  terms,  a  plasma  confined  by  a  magnetic  field 
must  be  in  a  state  of  non-equilibrium.  This  is  because  the  only  equi¬ 
librium  distribution  is  a  Maxwellian,  which  is  a  function  only  of  the 
Hamiltonian.  Thus  instabilities  produced  within  a  plasma  are  really 
the  result  of  the  plasma  trying  to  approach  an  equilibrium  configuration 
independent  of  any  magnetic  fields.  One  mechanism  for  attaining  equi¬ 
librium  is  through  binary  collisions.  However  there  are  mechanisms  by 
which  collisionless  plasmas  (collision  frequency  is  small  compared  to 
wave  frequency)  can  approach  equilibrium.  It  is  this  latter  type  of 
instability  which  will  be  of  concern  here. 

Instabilities  can  be  divided  into  two  types:  hydromagnetic 
instabilities  and  microinstabilities.  Hydromagnetic  instabilities  are 
those  derived  from  the  fluid  equations  where  essentially  the  magnetic 
field  lines  are  "frozen"  into  the  plasma  as  a  result  of  the  constancy  of 
flux  through  a  material  surface  moving  with  the  fluid.  These 
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instabilities  are  large  scale,  rapid  and  chaotic,  moving  the  plasma  as 
a  whole  away  from  a  confined  configuration.  The  microinstabilities 
however  can  only  be  derived  from  microscopic  equations  such  as  the 
Boltzmann  equation.  These  instabilities  develop  from  small 
perturbations  within  the  plasma  which  grow  with  time.  The  resulting 
turbulence  and  diffusion  ultimately  destroy  the  confinement  of  the 
plasma.  Only  small  groups  of  particles  may  be  initially  involved  with 
microinstabilities  as  compared  to  the  whole  plasma  in  the  hydro- 
magnetic  case. 

One  of  the  important  microinstabilities  which  has  been  investigated 

1  2  3 

is  the  collisionless  drift  wave  (Krall  ,  Rosenbluth  ,  Politzer  )  which 
arises  as  a  result  of  the  density  gradients  inherent  in  any  confined 
plasma.  The  wave,  which  is  driven  by  a  resonant  particle  interaction 
is  shown  to  be  characterized  by  a  low  frequency  (1-10  KHz)  and  a  small 
growth  rate.  It  was  found  that  in  many  plasma  devices  non-uniform 
electric  fields  were  present  which  resulted  in  a  non-uniform  rotation 
of  the  plasma.  Deviations  from  a  gaussian  density  profile  also 

4 

resulted  in  non-uniform  rotations.  The  paper  by  Capjack  and  Stringer 
investigated  the  effects  of  non-uniform  rotations  on  the  drift  wave 
when  the  density  was  sufficiently  high  so  that  terms  involving  the 
ratio  of  the  Alfve'h  velocity  to  the  velocity  of  light  could  be  ignored. 
Numerical  studies  were  made  showing  the  dependence  of  the  growth  rate 
and  frequency  of  the  drift  wave  on  the  magnitudes  and  scale  lengths 
of  the  electric  fields  and  density  gradients  present. 

The  purpose  of  this  thesis  will  be  first  to  provide  some  back- 
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ground  on  the  drift  mode  and  the  experimental  conditions  under 
which  it  can  be  observed.  A  brief  explanation  for  the  steady  state 
electric  fields  found  in  some  plasma  devices  will  be  given.  The 
drift  mode  will  then  be  examined  in  a  non-unif ormly  rotating 
cylindrical  plasma.  Terms  involving  the  ratio  of  the  Alfven  velocity 
to  the  velocity  of  light  will  be  retained  in  the  differential 
equation  describing  the  drift  mode.  This  will  allow  lower  density 
plasmas  to  be  investigated.  Approximate  analytic  solutions  both 
qualitative  and  quantitative  in  nature  will  be  given  using  the  WKB 
method  and  a  Weber-Hermite  approximation.  The  analytical  results 
will  help  clarify  the  numerical  results  presented  and  will  illustrate 
the  physical  mechanisms  involved.  Conclusions  will  be  drawn  concern¬ 
ing  the  effects  of  non-uniform  electric  fields  and  density  gradients 
on  the  drift  instability  for  density  regimes  typical  of  many  plasma 


devices . 
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CHAPTER  II 

THE  DRIFT  INSTABILITY 

2.1  Physical  Mechanism 

5  6  7 

Following  the  methods  of  Davidson  and  Kammash,  Stringer,  Meade 

O 

and  Sagdeev,  a  heuristic  derivation  of  the  frequency  and  growth  rate 
of  the  drift  wave  will  be  given. 

A  plasma,  if  it  is  to  be  in  equilibrium  according  to  the  fluid 
model,  must  exhibit  diamagnetism.  From  the  Navier-Stokes  equation  for 
ions  or  electrons  in  equilibrium  one  has: 

3  X  B  =  VP  =  VnkT  (2.1) 

where  J  =  current  density  =  neV 
B  =  magnetic  field 
n  =  density 
T  =  temperature 
P  =  pressure 

V  =  macroscopic  velocity 
e  =  charge 

k  =  Boltzmann’s  constant 

For  a  cylindrical  geometry  where  the  density  gradient  is  in  the  radial 
direction  and  a  uniform  magnetic  field  is  in  the  axial  direction, 
equation  (2.1)  yields: 
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where 


-kT  ne(r)'gQ 
eB  ne(r) 


-  -kT  n1 (r)  eQ 


eB  n  (r) 


electron  diamagnetic  drift  velocity  in  the 

azimuthal  direction  (e  ) 

0 


(2.2) 


ion  diamagnetic  drift  velocity  in  the 
azimuthal  direction 


♦ 

denotes  differentiation  with  respect  to  r 

i,e  denotes  ions  or  electrons  respectively 

The  electrons  and  ions  thus  have  diamagnetic  drift  velocities  given 
by  equation  (2.2)  which  provide  currents  necessary  for  hydrodynamic 
equilibrium.  By  using  a  simple  model  requiring  charge  neutrality 
when  ions  move  across  the  field  lines  with  an  electric  drift  velocity 
and  electrons  adjust  their  density  according  to  a  Maxwellian,  it  can 
be  shown  that  a  wave  may  exist  propagating  in  the  electron  diamagnetic 
drift  direction  with  a  phase  velocity  (V  )  equal  to  that  of  the 
electron  diamagnetic  drift  velocity.  Ion  motion  in  the  axial  direction 
is  neglected  since  it  will  be  shown  later  that  a  drift  wave  will  locate 
itself  in  a  region  where  the  ion  thermal  velocity  is  much  less  than  the 
phase  velocity  of  the  wave. 

One  begins  by  considering  an  inhomogeneous  neutral  plasma  in  a 
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uniform  axial  magnetic  field  with  a  density  gradient  in  the  radial 
direction.  An  electrostatic  wave  perturbation  (normal  mode)  is 
assumed  to  develop  with  a  wave  vector  K  in  the  axial  and  azimuthal 
direction  (Figure  1) . 

e 


Figure  1  -  Simple  plasma  model  in  a  local  co-ordinate  system 

The  perturbed  potential  is  of  the  form: 

<J)  =  exp  !  i(mO  +  K,,z  -  Wt)j  (2.3) 

where  K,,  and  m  are  the  wave  numbers  in  the  ^  and  eQ  directions 
respectively.  The  electric  field  corresponding  to  is  given  by: 

E  =  -V^  -K „<J>  ez  -  imc)>  eQ 


r 


(2.4) 
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The  equation  of  continuity  for  ions  is: 


Sn1  +  V» (n^V^)  =  0 

3t 


(2.5) 


which  is  linearized  by  letting 


i  i  i 

n  =  n  +  n. 

o  1 


V1"  =  V1  +  (since  the  ions  are  cold) 

oil 


where  is  the  first  order  E  X  B  drift  velocity  in  the  radial 

direction  resulting  from  the  perturbation  The  continuity  equation 
then  gives  an  expression  for  the  perturbed  ion  density  n^: 


n. 


,  l 

-m  o-  n 
1  o 

r  3W 


(2.6) 


With  the  flow  of  electrons  not  limited  in  any  way,  they  are  locally 
Maxwellian  along  the  perturbation  and  so  the  perturbed  electron 
density  is: 


0  0 

n,  =  ne  e<t>1 

1  o  J 

kT 


(2.7) 
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By  using  the  quasi-neutrality  condition  (n  -  n  ),  which  is  valid  for 
sufficiently  high  density,  it  can  be  established  that  the  dispersion 


relation  is: 


W  = 


-mkT 

eB 


n 


e 

n  r 

o 


=  mW. 


(2.8) 
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-kT 

e.B 


e 

n  r 

o 


the  electron  diamagnetic  angular 

0  I 

drift  velocity  (V^  j r) 


From  this  simple  model,  the  drift  wave  is  seen  to  be  neutrally 
stable  (imaginary  W  =  0)  having  an  angular  frequency  equal  to  a 
multiple  m  times  the  electron  diamagnetic  angular  frequency  (W  ) . 

hiN 

Since  this  model  neglects  microscopic  effects  by  ignoring  the  finite 
velocity  spread  of  the  particles,  it  cannot  correctly  predict  the 
growth  rate  of  the  wave.  There  are  mechanisms  limiting  the  flow  of 
neutralizing  current  along  the  field  lines  resulting  in  an  exchange 
of  energy  between  the  wave  and  the  particles.  These  mechanisms  take 
the  form  of  collisional  resistance  for  a  collisional  plasma  and 
resonant  particle  effects  for  a  collisionless  plasma.  The  result  is  a 
positive  imaginary  part  of  W  giving  rise  to  a  growing  unstable  wave. 
The  following  argument  (Meade^)  gives  an  estimate  for  the  growth  rate 
in  a  collisionless  plasma. 

A  wave  will  grow  in  a  plasma  when  particles  transfer  energy  to  the 

wave.  One  can  look  at  the  particle  wave  interaction  as  being  composed 

of  two  types  of  interactions.  One  interaction  results  from  particles 

travelling  slightly  faster  or  slower  than  the  phase  velocity  of  the 

9 

wave.  This  gives  rise  to  Landau  damping  (Stix  ).  The  other  inter¬ 


action  results  from  particles  travelling  exactly  at  the  phase  velocity 
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of  the  wave.  In  this  case  particles  can  either  give  energy  to  the 
wave  or  take  energy  from  the  wave,  depending  on  the  direction  of  wave 
propagation  and  sign  of  the  charged  particles.  The  physical 
mechanism  for  this  second  type  of  interaction  is  now  given. 


Figure  2 


Equi-potential  surfaces  of  a  collisionless 
drift  wave  in  a  local  co-ordinate  system  in 
cylindrical  geometry 
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Figure  2  shows  a  drift  wave  whose  phase  velocity  in  the 

azimuthal  direction  is  the  same  as  the  electron  diamagnetic  drift 

velocity.  The  streaming  ions  or  electrons  along  the  magnetic  field 

lines  interact  with  the  z-component  of  E(EM).  For  the  moment,  electrons 

only  are  considered.  At  position  "a"  resonant  electrons  experience  a 

constant  E  X  B  drift  in  the  radial  direction.  At  position  "b"  resonant 

—y  ~y 

electrons  experience  the  E  X  B  drift  in  the  opposite  direction.  The 
direction  in  which  resonant  electrons  travel  is  such  that  the  electron 
current  at  "a"  opposes  E,,  giving  energy  to  the  wave  while  the  current 
at  "b"  is  in  the  same  direction  as  E,,  taking  energy  from  the  wave. 
Because  of  the  density  gradient,  there  will  be  more  electrons  at  "a" 
than  at  "b"  and  thus  the  net  result  is  that  resonant  electrons  give 
energy  to  the  wave.  Since  the  ions  are  of  opposite  charge,  they  will 
take  energy  from  the  wave.  The  quantitative  aspects  of  the  above  are 
now  given. 

The  net  power  per  cycle  delivered  to  the  wave  by  the  resonant 
electrons  and  ions  is  shown  to  be^: 


(2.9) 


' 
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where  the  assumptions  made  are  that: 

i.  The  velocities  of  the  ions  and  electrons  in  the  axial 
direction  follow  Maxwellian  distributions  which  are  at 
the  same  temperature  T. 

ii.  The  phase  velocity  of  the  wave  in  the  axial  direction 


(W/K,,)  is  much  less  than  the  ion  thermal  velocity  (V^  ) 

I ri 


0 

or  the  electron  thermal  velocity  (V  ) 

i  H 


iii.  The  electron  and  ion  densities  are  approximately 
equal  (n1  -  n6  =  n) 


The  net  power  per  cycle  taken  from  the  wave  by  electron  and  ion 

9 

Landau  damping  is  given  by  : 


(2.10) 


where  the  same  assumptions  are  made  as  for  (2.9). 


The  energy  per  unit  volume  in  the  drift  wave  is  given  by: 


(e  =  dielectric  tensor  (2.11) 


where  the  assumptions  made  in  addition  to  those  of  equation  (2.9)  are 
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i.  the  wave  frequency  W  is  much  less  than  the  ion  cyclotron 
frequency  (fh) 

ii.  the  ion  Larmor  radius  (a^)  is  much  less  than  r/m 

Since  E  varies  as  exp(-iWt)  =  exp (iW t)exp(W  t)  where  W  is  the 

K  1  R 

real  part  of  W  and  W  is  the  imaginary  part  of  W,  the  equation  for 
the  conservation  of  energy  is: 


2W  (e .E)*E 
8tt 


(J-E)r  + 


(2.12) 


When  equations  (2.9),  (2.10)  and  (2.11)  are  substituted  into  (2.12) 


the  expression  for  the  growth  rate  W^.  becomes: 


+  W)  ./'  W  ^ 


lK"VTH 


exp 


(2.13) 


It  can  be  seen  from  (2.13)  that  if  ion  Landau  damping  is  allowed 

(W/K,,V^tt  <<  1)  the  drift  mode  with  W  -  mW  will  be  highly  damped 

1  rl  K  EN 

(since  v1  <<  Ve  )•  If  the  drift  mode  is  in  a  region  where  ion 
TH  TH' 

Landau  damping  does  not  occur  (W/K„V  H  >>  1)  the  growth  rate  becomes 


* 


, 


(2.14) 


One  can  see  that  in  this  case  there  is  a  competition  between 
the  process  of  electron  Landau  damping,  which  causes  the  wave  to  damp, 
and  resonant  particle  interaction  which  causes  the  wave  to  grow.  By 
using  the  fluid  model  it  was  found  that  W  =  mW  ;  it  might  be  con¬ 
cluded  that  the  wave  is  neutrally  stable  since  the  two  interaction 
processes  balance  each  other.  It  is  thus  evident  that  any  effect 
which  causes  W  <  W  will  give  rise  to  an  instability  in  the  region 
where  ion  Landau  damping  can  be  ignored.  In  chapter  three  it  will  be 

shown  that  the  finite  Larmor  radius  of  the  ions  allows  W  <  W  . 

rAN 

One  can  interpret  the  physical  mechanism  of  the  drift  mode  in  a 

6  8 

slightly  different  manner  (Stringer,  Sagdeev  )  which  shows  the  two 
types  of  interactions  mentioned  above  really  to  be  manifestations  of 
the  same  basic  Landau  interaction.  The  Landau  theory  shows  that  if  the 
number  of  particles  moving  slightly  slower  than  the  phase  velocity  of 
the  wave  exceeds  the  number  moving  slightly  faster,  the  plasma  will  be 
stable  against  self-excitation.  For  a  homogeneous  plasma  this  condition 
becomes  3f  /9V  TT  T7z  <  0  which  is  satisfied  for  a  Maxwellian 

0  2  VVPH 

2 

distribution.  (f  is  the  equilibrium  distribution  function,  V  is 

0  r  H 

the  phase  velocity  in  the  axial  direction  and  is  the  particle 
velocity  in  the  axial  direction.) 


* 
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However,  in  an  inhomogeneous  plasma,  terms  involving  Vf  occur  in 

o 

the  Vlasov  equation  (Appendix  A) .  The  stability  criteria  then 

becomes  [ 3f/3V  -  (m/r) (V  /Wft  )3f  /3r]!  Trz  <  0.  The  term 

L  o  z  z  e  o  J  1 V  =Vn„ 

z  PH 

-(m/r)(V  /Wfi  )3f  /3r  i-s  positive  for  a  wave  propagating  in  the 

Z  G  0 

electron  diamagnetic  drift  direction.  If  the  plasma  has  a  Max¬ 
wellian  distribution  where  the  thermal  velocity  of  the  particles  is 
much  greater  than  the  phase  velocity  of  the  wave,  the  term 

will  be  a  small  negative  quantity  and  thus  the  stability 
criteria  can  be  violated.  What  has  happened  is  that  the  density 
gradient  has  upset  the  number  of  particles  going  slightly  slower 
or  faster  than  the  wave.  Resonant  particles  are  brought  in  from 
regions  of  different  densities  distorting  the  velocity  distribution 
The  number  of  particles  going  slightly  faster  than  the  wave  then 
exceeds  the  number  going  slightly  slower  giving  rise  to  instability 
The  following  conclusions  can  now  be  drawn  concerning  the  drift 
mode : 


i.  The  drift  wave  occurs  in  a  frequency  range  such  that 
Ti 


V"  <<  vz  <<  ve  . 

TH  PH  TH 


The  inequality  >>  V^u  is 

r  n  1 


necessary  to  avoid  ion  Landau  damping  which  will  tend 

Z  G 

to  stabilize  the  wave.  The  inequality  << 

is  necessary  so  that  the  density  gradient  can  distort 
the  distribution  function  sufficiently  to  allow  the 
number  of  particles  moving  faster  than  the  wave  to 


exceed  the  number  moving  slower. 


* 


' 
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ii.  Any  effect  which  allows  the  wave  frequency  to  be 
less  than  the  electron  diamagnetic  frequency  will 
give  rise  to  instability. 

iii.  The  growth  rate  W^.  is  much  less  than  the  frequency 

W  (equations  (2.8),  (2.13) 

K 


2.2  Experimental  Conditions 

3  10 

Experimental  identifications  of  the  drift  mode  5  have  been 
successful  in  devices  known  as  "Q-machines" .  A  typical  "Q-machine" 
consists  of  a  cylindrical  vacuum  approximately  100  centimetres 
long  and  5  centimetres  in  diameter  in  which  a  uniform  axial  magnetic 
field  of  about  2000  to  5000  Gauss  is  produced.  At  each  end  of  the 

cylinder  are  tungsten  hotplates  heated  to  temperatures  of  about 
2000  degrees  Kelvin  by  electron  bombardment  from  a  hot  filament. 

One  plate  is  usually  free  to  move  so  that  the  axial  length  of  the 
machine  can  be  varied.  The  plasma  is  produced  thermionically  by 
spraying  a  beam  of  neutral  lithium  or  sodium  atoms  into  the  hot¬ 
plates.  The  beam  then  becomes  ionized.  Thermoelectric  emission  of 
electrons  from  the  hot  plates  provide  the  charge  necessary  to  keep 
the  plasma  "quasi-neutral".  Due  to  the  random  nature  in  the  method 
of  plasma  production, it  is  assumed  that  the  electron  and  ion 
velocity  distributions  in  the  axial  direction  are  Maxwellians 


■ 

- 
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having  a  temperature  equal  to  that  of  the  hotplate.  Typical 

8  3 

particle  densities  are  of  the  order  10  /cm  which  puts  the  plasma 
in  a  collisionless  regime  since  the  mean  free  path  is  several 
machine  lengths. 

Two  important  aspects  in  "Q-machines"  are  the  density 
gradients  and  radial  electric  fields  present.  Density  gradients 
must  be  present  when  a  plasma  is  confined.  A  typical  variation 
of  ion  density  with  radius  is  shown  in  figure  3a.  The  fact  that 
the  ion  density  profile  does  not  exactly  fit  a  gaussian  gives 
rise  to  non-uniform  ion  and  electron  angular  diamagnetic  velocities 
(figure  3b) .  Time  independent  radial  electric  fields  can  arise 
for  several  reasons.  Since  there  is  no  accumulation  of  charge  at 
the  end  plates,  the  thermal  efflux  of  electrons  from  the  plasma 
into  the  end  plate  sheath  must  be  equal  to  the  thermoelectric 
emission  of  electrons  from  the  end  plate.  The  flux  of  electrons 
from  the  plasma  to  the  plate  is  given  by: 


J 


F 


1/4  n6 

o 


(2.15) 


The  electron  flux  emitted  by  the  plate  is  given  by  the  Richardson 
current : 


exp 


(2.16) 


where 


the  work  function  for  the  end  plate 
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n1xlO 

o 

,  -3x 

(cm  ) 


1 

rn 

o 


. 
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ip  =  the  potential  of  the  plasma  with  respect  to  the  end  plate 


1.0  x  10  electrons  cm 


A 


By  equating  J  with  J  ,  solving  for  <f>  and  taking  its  gradient,  one 

Ij  r 


obtains^: 


(2.17) 


Thus  density  gradients  in  the  plasma  as  well  as  temperature  gradients 
in  the  end  plates  give  rise  to  electric  fields.  Since  the  gradients 
are  predominantly  radial,  the  electric  field  will  be  radial.  The 
diffusion  and  mobility  of  the  ions  and  electrons  can  give  rise  to 
particle  losses  which  result  in  electric  fields.  The  electron  and 
ion  particle  currents  are  given  by: 


T  e  ->■ 

e  -  -nU  «  E 


r. 

l 


(2.18) 


where 


electron  (ion)  mobility  tensor 


electron  (ion)  diffusion  tensor 


e  i 


n  (particle  density) 


n  -n 


By  putting  f 


1.  (steady  state)  one  obtains: 


e 


■ 
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-> 


E 


+ 


u1)"1 


D 


(2.19) 


Since  electrons  diffuse  more  rapidly  than  ions,  an  electric  field 

will  develop  (equation  2.19)  which  will  slow  down  the  electrons, 
accelerate  the  ions  and  thus  allow  the  two  particle  currents  to 
become  equal.  Also,  radial  electric  fields  are  often  imposed 
externally  on  the  plasma  for  purposes  of  investigation.^  The 
end  plate  is  split  into  two  concentric  parts  separated  by  a 
small  gap.  By  applying  a  voltage  between  these  two  sections, 
an  adjustable  radial  electric  field  can  be  obtained.  A  typical 
variation  of  voltage  with  radius  is  given  in  figure  4a  with  the 
resulting  angular  electric  drift  velocity  profile  (figure  4b) . 

In  the  type  of  "Q-machine"  just  described  one  observes 

3 

localized  oscillations  near  the  region  of  maximum  density  gradient. 
These  oscillations  propagate  with  an  azimuthal  velocity  of  the  same 
order  of  magnitude  as  the  electron  diamagnetic  drift  velocity.  A 
standing  wave  exists  in  the  radial  direction  (usually  a  zero’ th 
radial  mode)  and  in  the  axial  direction  (  parallel  wave  length  is 
approximately  twice  the  machine  length) .  The  radial  electric  fields 
show  a  stabilizing  influence  when  the  electric  rotation  opposes  the 


electron  diamagnetic  rotation 


> 


, 

. 
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Figure  4a  -  Variation  of  voltage  with  radius 


Figure  4b  -  Variation  of  1  with  radius 

r  dr 
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2.3  Theoretical  Model 


Any  theoretical  efforts  to  explain  the  properties  of  the 

drift  mode  must  be  modeled  to  the  experimental  situation 

described  in  the  above  section.  This  is  accomplished  by  choosing 

suitable  functions  to  represent  the  profiles  shown  in  figures  3 

4 

and  4  (Capjack  and  Stringer  ).  The  profile  of  figure  3b  is 
represented  as: 


I 


( 1  +  bexp 


(r-rM)2/LN 


) 


L  r  (1  +  b) 
N  M 


(2.20) 


where : 

L  =  the  density  scale  length 
n 

r  ,  =  the  radius  at  which  a  is  maximum 

M 

b  =  the  parameter  representing  the  non-uniformity  in  a 
The  resulting  ion  density  profile  (figure  3a)  is  given  by: 


l 

n 

o 


n  exp 
o 


-1 


r  +  bL 


-  exp 


VM(1+b) 


(r~rM)2/LN 


N 


LN  exp 


[  (W  ] 


(r-rM)/LN 

r  2 

+  r^  |  exp(-z')  dz 
-rM/LN 


(2.21) 
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Finally,  the  profile  of  figure  4b  is  represented  as: 


II 

"it 

o 

10  Max 

I  it 

lr  7 

r* 

r  3r 

volts /cm 


2 


sgn 


(2.22) 


where  L  =  the  electric  scale  length 

Cj 

It  will  be  assumed  that  a  and  (1/r  3 / 3 r )  are  maximum  at  the  same 
radius  r^.  The  electron  diamagnetic  frequency  (W^)  ion  dia¬ 
magnetic  frequency  (W„) ,  electric  drift  frequency  (  W  )  and  the 

Ej 

0 

electron  density  (n^)  can  all  be  determined  from  equations  (2.20), 
(  2.21)  and  (2.22).  These  are  given  as  follows: 


(2.23) 


W 


kT  a  = 
eB 


"Vl+  beXP!"(r'rM)2/LN 
(1+b) 


) 


(2.24) 


e 

n 

o 


n 


1  +  (CA/C)' 


(  2We  +  rWE,)/fiij 


(See  A. 4) 


(2.25) 


W. 


EN 


=  -W. 


/kl\ 

M2 

*(2We  +  WE'x-ar)  +  (3WE*  +  W£  r) 

7 

U±  +  (Ca/C)2(2We  +  WE'r) 

(2.26) 
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where 

CE  =  Max  (WE) 

CNE  =  Max(Wj 

CA  =  Alfven  velocity  (Appendix  A,  equation  A. A) 

All  the  properties  of  the  unperturbed  plasma  are  assumed  to  be 
adequately  described  by  equations  (  2.20  -  2.26).  Table  1  gives 

the  typical  numerical  values  for  the  important  plasma  parameters 
found  in  "Q-machines " .  All  subsequent  approximations  made  are 

based  on  these  values. 
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Table  1  -  Plasma  Parameters 

electron  density . 

ion  density . 

ion  Larmor  radius . 

electron  Larmor  radius . 

electron  or  ion  temperature . 

magnetic  field  (uniform) . 

peak  diamagnetic  rotation . 

density  scale  length . 

non-uniform  rotation  parameter. . . 

radius  at  peak  rotation . 

electron  cylotron  frequency . 

ion  cyclotron  frequency . 

parallel  wave  number . 

ion  thermal  velocity . 

electron  thermal  velocity . 

plasma  radius . 

peak  electric  rotation . 

electric  scale  length . 

electron  mass . 

ion  mass  (Lithium  plasma) . 


in8/  3 
10  /cm 

in8/  3 
10  /cm 

5x10  ^  cm 
-4 

4x10  cm 

2000  °K 
3000  Gauss 

4 

2x10  /sec 

4x10  ^cm 

10 

8x10  ^cm 

5x10  "^/sec 

4xl0^/sec 

_? 

3x10  /cm 
2xl0^cm/sec 
2xl0^cm/sec 
2.5  cm 


~  C 


NE 


m  L. 


N 


9.1x10  ^gm 


11 . 7x10 


-24 


gm 


■ 


, 
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CHAPTER  III 


ANALYSIS  OF  THE  DRIFT  MODE 


3.1  Differential  Equation 

The  differential  equation  describing  the  drift  instability 
in  the  presence  of  non-uniform  radial  electric  fields  is 
derived  in  Appendix  A  for  the  variable  ip ,  which  is  the  perturbed 
guiding  centre  displacement.  The  equation  is: 


(A.  1  ) 


By  making  the  substitution  \p  —  S^y  in  the  above  equation,  one 


obtains : 


y 


M 


1  Q(W,r)y  =  0 

2 


(3.1) 


a . 
i 


where : 


y(r) 


3  i  9  ' 

w  r  n‘  (W  -  mWA  +  W  (CA/C)  ) 
o  o  o  o  A 


~'2 


\  °  /  L 

2(w  -mW  )  Z  (W  /K„V®  )  +  2  (W  -mWjZ(W  /K,,V^  ) 
oENoIH  o  olH 


+  a^FL(W,r) 


W  (1-K  CA/C)2-mW,v/W  ) 
o  A  o 


Q(W,r)  = 


. 
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m2-l  (CA/C)  2K2-(W/W  )2a  1  /S  ’\  '  1  fs 

„  .  A  o  o  o 

W  » r )  “  o  +  o  +  i  +  I 


r2  (  l+(C./C)2-mWA/W  )  2  Is 

A  o  \  o, 


4  \S 


o  I 


2  2 

One  can  note  that  a.FT  is  of  the  order  (a./L)  times  that  of 

l  L  l 

the  first  or  dominant  term  where  L  is  a  characteristic  length 

(density  or  electric  field  scale  length)  in  the  plasma.  The 

2 

term  a.FT  is  thus  considered  as  a  correction  term  associated 

l  L 

with  finite  Larmor  radius  effects.  The  remainder  of  this  chapter 

will  deal  with  the  solution  of  (3.1)  for  the  real  parts  (W  ) 

K 


and  imaginary  parts  (W^)  of  the  complex  frequency  W. 


3.2  WKB  Solution,  Local  Approximation 

It  is  assumed  that  Q  has  only  two  turning  points  z  and  z ^ 

where  Q(z  )  =  Q(z?)  =0.  In  order  that  a  bounded  solution  exists 
with  the  solution  far  from  the  turning  points  "joining  smoothly" 
with  the  solution  in  the  region  of  the  turning  points,  the  follow 
ing  quantization  condition  must  hold  in  the  complex  plane: 


1 

a . 
i 


dz 


(n  •+  %)tt  (n  is  an  integer) 


(3.2) 
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If  the  solutions  y(W,z)  are  to  be  analytically  continued  to 
the  real  "r"  axis  such  that  y(W,r)  0  as  r  ->  00 ,  the  branch  cuts 
from  z_  and  z„  must  join  at  z  =  00  and  the  Anti-Stokes  lines 
must  cross  the  real  axis  (figure  5).  The  latter  condition  will 
always  hold  since  Q  is  almost  real  with  almost  real  turning 
points . 


Im  z 


Figure  5  -  Branch  cuts  and  Anti-Stokes  lines  characteristic 
of  a  valid  WKB  solution.^ 

One  now  assumes  that  highly  localized  modes  exist,  implying 
that  the  turning  points  will  be  close  together.  Thus  if  Q  is 
expanded  in  a  Taylor  series  about  the  point  z^  where  Q'(zq)  =  0, 
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keeping  terms  to  order  (  z-z  )  only  will  result  in  a  reasonable 

o 

approximation  for  Q  in  the  region  between  the  turning  points. 

If  the  expansion  for  Q  is  substituted  in  equation  (3.2)  the 
result  is: 

Q(z  )  =  3  (n  +  %)tt  (-Q"a2)  2 

°  8  1  (3.3) 

The  local  approximation  (Krall^")  involves  setting  the  right 

hand  side  of  equation  (3.3)  equal  to  zero  on  the  assumption  that 
2  t  2 

(-Q"a.)2  ^  0(a./L^T)  <<  1.  The  dispersion  relation  then  becomes: 

Q(zq,W)  =  0  ,  Q’(zo,W)  =  0  (  3.4) 


The  function  Q(zq,W)  can  be  written  as  f(z^,W)  +  ig(zQ,W).  If 

the  assumptions  that  |g|  <<  If  I  and  W  <<  W  are  used,  equation  (3.4) 

IK 

can  be  written  as: 

\ 

Q(zq,W)  -  f  (WR)  +  ig(WR)  +  i3f  WI  =  0  (  3>5) 

3W 

The  dispersion  relation  then  becomes: 


f(V 


o 

1  Z=Z 

o 


-g  (WR) 

~TT 

3W  w=wn 

K 


Z  =  Z 

o 


(  3.6) 


■ 
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Based  on  the  physical  arguments  of  chapter  two,  equation 

(3.1)  will  only  be  considered  in  the  limits  V^TT  <<  W  /K..  << 

TH  o  i 

13 

The  Z  functions  can  then  be  written  as: 


z(wo/k„v®,h)  =  1  +  iAwo/K„v®H  (Wo/K"V1H)  <<:  x 

z(uo/k„vjh)  =  -%(v^hk„/wo)* 2  (  wo/k„v^h)  »  1 


Equation  (3.6)  is  now  used  with: 


2(Wo-mWEN)  -  (vjHK„/Wo)  (W^nW^)  +  a.F^ 


f  = 


l+(C./C)2-mW,/W  1 
A  o 


W 

o 


1+(CA/C)  -mW./W 
A  o 


2/TT(Wo/K„VTH)(Wo-mWEN) 


g  = 


w 


o 


1+(C  /C)  -mW./W 

A  O 


2  ^  2 

By  keeping  terms  to  first  order  in  a^^  and  (V,rHKM/W  )  and  assuming 
quasi-neutrality , one  obtains  the  following  expressions  for  the 
frequency  and  growth  rate  of  the  drift  wave: 
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WR  =  ”% 


2  * 

1  -  aTFT  + 
l  L 


(vthk"/wr>2 


r=r. 


M 


wi  ■  (H/K"vra)(mWEN-WR> 

“  /,r(nWEN)2  aiFL  -  <VTHK"/WR>2  ! i — 

K„Ve 

"  TH 

(3.7) 

where 


W  -  mW„ I 
E 


+ 


1  WEN 

4  V 


WE  z'3+2  (cA/e)2 

v(  4 


+ 


2  / 1+(CA/C)2  \ 
\uh(ck/c)2} 


The  result  for  the  growth  rate  W  Is  basically  the  same  as 

that  obtained  in  chapter  two  (equation  2.14)  except  the  wave 

frequency  TJ  is  replaced  by  W°  the  wave  frequency  in  a 
R  R 

reference  frame  rotating  at  the  peak  electric  angular  frequency  C  . 

E 
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The  effect  of  the  radial  electric  field  is  to  "doppler-shif t" 

the  wave  frequency  W  .  The  finite  Larmor  radius  effects,  represented 

R 

2  * 

by  the  terra  a.F  give  rise  to  a  positive  growth  rate.  Ion  Landau 
1  .Li 

damping  is  seen  to  have  a  stabilizing  influence  through  the  term 


,oN  2 


(  V  K„/(L)  •  Since  the  dominant  term  (m  ~h) /r  in  F  is  magnified 
1  li  K  1 


/  2 


2  . 


as  (C  /C)  increases,  a  lower  density  plasma  should  exhibit  a 
higher  growth  rate.  It  will  be  seen  in  the  next  section  that  the 
local  approximation  cannot  adequately  describe  the  effects  of  a 


non-uniformity  in  the  electric  and  diamagnetic  rotations.  This  is 
2  L 

because  (Q"a~) 2  turns  out  to  be  of  the  order  (a_^/L^)  rather  than 
2 

(  a./L)  when  non-uniform  effects  become  important.  The  local 
l  N 

approximation  then  becomes  invalid. 


3.3  Weber-Hermite  Approximation 

Equation  (3.1)  will  be  approximated  by  a  Weber-Hermite 

equation.^  It  will  be  shown  that  because  of  non-uniform  rotations, 

the  finite  larmor  radius  affects  the  growth  rate  to  order  (eu/L  ) 

2 

rather  than  .(^/L^  as  in  the  local  approximation. 

In  the  region  V*  «  W  /K„  <<  V®R,  Q  (W,r)  can  be  written  as: 


Q(W,r) 


2«o  -  mWEN)  (1  +  i/lTWo/K"VTH) 

V1  +  (CA/C)2  +  mWEN/Wo? 


+  a^FL(W,r) 


(3.8) 


, 
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where  the  ion  term  has  been  ignored  since  it  only  contributes 

a  very  small  stabilizing  influence  in  the  region  of  interest. 

■k 

Q  is  expanded  about  the  values  r=r  and  W=W  ^mC^^+mC^  where  terms 

M  NE  E 

3  *  3 

of  order  (r-r  )  ,  (W-W  )  and  higher  are  ignored. 


k  0  k  0  0  ( V  ■ 

Q(W,r)  =  Q(W  ,r  )  +  3Q(W-w")  +  3  Q(W-W  )  +  3_Q  _ _ M; 

»  1?  2  3r2  2 


(3.9) 


The  assumptions  made  in  equation  (3.9)  are: 

i.  The  values  of  W°  to  be  obtained  are  close  to  imC  , 

R  ME 

the  peak  electron  diamagnetic  frequency.  Thus  the 

terms  (W°-mC.T„)n  Qn(W)/n!  can  be  ignored  for  n  >  3. 

K  Mb 

Also,  the  growth  rate  W  is  assumed  to  be  much  less 

than  W°, 

K 

li .  The  modes  are  localized.  Thus  the  magnitude  of  the 

eigenfunctions  are  appreciable  only  in  the  region 

(r-r  ) /L  <<  1  where  the  expansion  in  equation  (3.9) 
M  N 

will  be  valid. 

iii .  Q  has  only  two  turning  points. 

iv.  Q  is  essentially  an  even  function  of  (r-r,,)  . 

M 
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v.  The  term  a.FT  (W,r)  is  made  constant  at  the  value 
l  L 

2  * 

a. F  (W,r).  Since  this  term  is  considered  a  small 

1  Li 

correction  to  Q,  its  variation  within  the  region 
of  localization  can  be  ignored  without  greatly 
affecting  the  final  results. 

By  evaluating  the  terms  in  equation  (3.9)  one  obtains: 


Q(W,r)  -  (1  +  i/mCNE/K„V®H) 


(w-w  ) 


(mCNE)(l+MCA/C)  ) 


+  aiFL(W’V 


(1  +  i/™cne/k"vto> 

1 

(mCNE)2(l+%(CA/C)2) 

__  -J 

2(1+%(Ca/C)2) 

1 


*  2 
(W-W  ) 


(l+i/TrmCNE/KnVaH) 


+ 


r 


(1  +  i/TrmCNE/KMVaH) 
(1-^(CA/C)2) 


b 


1+b 


+ 


W 


(3.10) 


The  differential  equation  (3.1)  can  now  be  written  as: 


+ 


2  2 
X  -  AZ(r-rM) 


0 


(  3.11) 


. 


where 


(1  +  i/™CNE/K„V®H)  (W*-W) 
ai  (mCNE)  (  1+'^CA/C2> 


+ 


(  1  +  i'/™CNE/K„VaH) 

-1 

1 

+ 

(mCNE)2(l+^(CA/C)2) 

2(1+^(Ca/C)2) 

(l+i/7rmCNE/KnVaH) 

2 
a . 
1 


A_  (1  +  1  m<WK"vra> 

(1+%(CA/C)2) 

Equation  (3.11)  is  just  the  Weber-Hermite  equation.  In  order 
that  y  -*  0  as  |r-r„|  ->  ±°°  (condition  for  localized  modes)  the 

i  M  i 

following  condition  must  be  satisfied. 

X  =  (2n+l) A  (3.12) 

When  the  expressions  for  X  and  A  are  substituted  in  equation  (3.12) 

e  2 

and  terms  of  order  (mC1.1_1/KIIV_1J)  are  ignored,  the  following 

Mb  1  n. 

expressions  for  the  "doppler-shif ted"  frequency  (W°)  and  the 
growth  rate  (W^)  are  obtained: 


b  +  ^CE^CNE^ 
,l+b  (LE/LN)2 


aiLNy 
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Figure  6  -  Growth  rate  and  frequency  versus  the 

important  plasma  parameters  represented 
through  the  variable  P.  (See  equation  3.13) 
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1  ~  v/l+E  (1-17/  2  D  ) ' 

2(  1-1/ 2D) 


- 

1 


V*w2 

(3/2-1/ 4D)P 

/l+P  (1-1/ 2D)1  -1) 

WvtJ  2 

2  (1-1/2D)  /l+P  (1-1/2D)' 

(1-1/ 2D) 2 

(3.13) 


where 


P  =  4(2n+l) (a./LN) 


(i+%(cA/c)  )/b_+(VV 


1+b 


(VV 


2  2  * 
+  4a.  Ft 
i  L 


D  =  (1-Ms(Ca/C)2) 


Hie  variable  P  contains  the  important  plasma  parameters  which 

determine  the  frequency  and  growth  rate  of  the  drift  wave.  Figure  6 

shows  W°  and  W  as  functions  of  P.  Only  the  shaded  region  of 

figure  6  will  be  of  interest  for  the  following  reasons: 

i.  The  approximations  made  in  equation  (3.10)  become 

invalid  as  P  becomes  of  the  order  unity. 

ii.  Experimentally,  the  modes  of  interest  are  for 

W°  -  mC  „  ,  which  implies  that  P  <<  1. 

R  NE 

In  the  limit  that  P  <<  1,  the  following  expressions  for  the  frequency, 
growth  rate,  eigenfunction  and  "potential  well"  Q  pertaining  to  the 


37 


drift  mode  are  obtained: 


(2n+l) 


(i-^s(cA/c)  )  b  + 

\l+b 


(C  /C  )\ 

V  E'  NEy 

(VV7_ 

(3.14) 


WI  ’  A(mCNE>2  , 

(2n+l)  (a±/LN) 

(l^(CA/C)2)(b  +  (ce/cne)^ 

U  2  * 

+  VlA 

2k"v?h 

\1+b  <VV7 

(mCNE 


WR> 


K„V, 


TH 


(3.15) 


yn(r)  =  "n  f/2A  (r_rM)]  exp["A(r"rM)2/^]  (3.16) 

where  H  are  Hermite  polynomials 


Q(W,r)  = 


•(l+i/TnnCNE/KMV®H) 


2.3s 


(1+^(CA/C)  ) 


(2n+l) 


/a.\  [  b  +  (CR/CNF) ,  2 


+  v  E  NEy 

v1+b  W 


+  a .  Ft  -  iW 
l  L  I 


+ 


(1  +  i/TrmCNE/KMVTH) 

(i+^(cA/c)2 


b  +  (VW 

x+b  (LE/LN)2 


(3.17) 
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The  important  feature  in  the  above  equations  is  that  non¬ 
uniformities  in  the  electric  and  diamagnetic  rotation  profiles, 

- — . 

represented  by  the  term  lb  (C_, /C_  )  have  a  far  greater  affect 

«  b  Nh 

1+b  (LE/LN)2 

on  the  growth  rate  of  the  wave  than  does  the  normal  destabilizing 

2  * 

term  a.FT  found  in  the  local  approximation.  This  is  because  the 
1  L 

2 

growth  rate,  which  depends  on  terms  of  the  order  (a^/L^)  when 
non-uniformities  are  absent,  now  depends  on  terms  of  the  order 
^ai^bP  (equation  3.15). 


3.4  Numerical  Solution 

The  numerical  solution  of  equation  (3.1)  was  obtained 
using  the  method  of  Appendix  B.  The  numerical  results  given  here 
illustrate  the  analytical  results  of  equations  (3.14  -  3.17)  for 
the  frequencies,  growth  rates,  eigenfunctions  and  "potential  wells" 
pertaining  to  the  drift  mode.  The  regions  for  which  Q  has  four 
turning  points  are  also  investigated.  Only  the  first  radial  mode 
will  be  of  interest  (n=0)  since  this  is  the  one  that  is  usually 
observed  experimentally.  All  parameters  are  given  as  in  Table  1 
unless  otherwise  specified. 


’ 
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Figure  7  -  The  doppler-shif ted  frequency  and  the  growth  rate  are  shown 

as  functions  of  (L  /L)  for  peak  ion  densities  of  5x10^ 

8  _3  E  N 

and  5x10  cm  .  (C„/C  )  takes  on  the  values  +.8,  -.8  and 

£j  Nt, 

0.  The  azimuthal  mode  number  m=l.  < - h depicts  the 

region  for  which  Q  has  four  turning  points. 
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Figure  8 


The  doppler-shif ted  frequency  and  the  growth  rate  are  shown 

as  functions  of  (C  /C  )  for  peak  ion  densities  of  5x10  ^ 

8  _  3  ^  ^ 

and  5x10  cm  .  (L  /L  )  =  .8.  The  azimuthal  mode  number 
m=l. 
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Figure  9 


-  The  doppler-shifted  frequency  and  growth  rate  are  shown 
as  functions  of  density  for  (CE/CNE)  =  ±-8>  (LE/LN)=*8‘ 
The  azimuthal  mode  number  m=l. 


VmCE 


mC 


NE 


W„ 


mC 


NE 


-1.40 


-0.84 


-0.28 


0.28 


0.84 


(LE/LN)Sgn(CE/CNE) 


1.1)0 


Figure  10  -  The  doppler-slrif  ted  frequency  and  the  growth  rate  are  shown 

as  functions  of  (L^/LO  for  peak  ion  densities  of  5x10  ^ 

and  5x10 ^  cm  ^ .  ( C  / C  )  takes  on  the  values  +.8,  -.8  and 

E  NE 


0.  The  azimuthal  mode  number  m=3 


->  depicts  the 


region  for  which  Q  has  four  turning  points. 
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Figure  11  -  The  doppler-shif ted  frequency  and  the  growth  rate  are  shown 

as  functions  of  (C  /C  )  for  peak  ion  densities  of  5x10 

h  uu 

and  5x10  cm.  (L  /L)  =  .8.  The  azimuthal  mode  number 
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Figure  12  -  The  doppler-shif ted  frequency  and  growth  rate  are  shown 

as  functions  of  density  for  (C  /C  )  =  ±.8,  (L  /L  )=.8. 
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The  azimuthal  mode  number  m=3. 
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Figure  13  -  The  real  parts  of  Q 
radius  r.  The  peak 


radius  (cm) 

and  y  (r)  are  shown  as  functions  of  the 

ion  density  n1  =  5xl08cm_3.  (C_/C.T_) 
o  E  NE 


=  .3,  (Le/Ln) 


=  .8 
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Figure  14  -  The  real  parts  of  Q  and  y  (r)  are  shown  as  functions  of  the 
radius  r.  The  peak  ion  density  n^  =  5x10 ^cm  ^ .  (C  /C  ) 


=  .8,  (Le/Ln)  =  .8 
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Figure  15  -  The  real  parts  of  Q  and  y  (r)  are  shown  as  functions  of  the 

i  o  —  3  ( C  /  C  } 

radius  r.  The  peak  ion  density  n  =  5x10  cm  E'  NE' 

J  o 


=  -.8,  aE/LN)  =  .8 


The  above  diagram  depicts  the  eigenfunction  and  the  potential  well  Q 
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when  the  growth  rate  is  minimized. 


n 

o 


48  - 


Figure  16  -  The  real  parts  of  Q  and  y q ( r )  are  shown  as  functions  of  the 

radius  r.  The  peak  ion  density  n^  =  5x10 ^cra  ^ .  (C  /C  ) 

o  E  NE 

=  -.8,  (Lg/Ljj)  =  .8 

The  above  diagram  depicts  the  eigenfunction  and  the  "potential  well"  Q 


when  the  growth  rate  is  minimized.  However,  the  plasma  density  is 
reduced  from  that  of  figure  15. 
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Figure  17 


The  real  parts  of  Q  and  y^(r)  are  shown  as  functions  of  the 

i  8  “  3 

radius  r.  The  peak  ion  density  n  =  5x10  cm  .  (C  /C,T10 

r  o  E  NE 

=  -.8,  (L/L,T)  =  .4.  It  can  be  noted  that  Q  has  four 
n  N 


turning  points 
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Figure  18  -  The  real  parts  of  Q  and  y  (r)  are  shown  as  functions  of  the 

_l  7  3 

radius  r.  The  peak  ion  density  n  =  5x10  cm  .  ( C  /C:  ) 

o  E  NE 


=  -.3,  (Le/Ln)  =  .A 


The  stabilizing  region  (region  1)  has  been  reduced  from  that  of 
figure  17  as  a  result  of  a  decrease  in  density. 
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3.5  Interpretation  and  Discussion 


From  the  previous  sections  it  is  seen  that  the  growth  rate 
depends  strongly  on  the  non-uniformity  of  the  diamagnetic  rotation 
and  on  the  manner  in  which  the  electric  rotation  interacts  with  the 
diamagnetic  rotation.  The  electric  rotation  has  little  effect  on 
the  growth  rate  when  the  electric  scale  length  is  much  greater  than 
the  density  scale  length.  When  the  electric  scale  length  becomes 
less  than  or  equal  to  the  density  scale  length,  the  growth  rate 
can  be  strongly  affected  by  the  electric  rotation.  This  is 
physically  reasonable  since  the  velocity  shear  resulting  from  the 
non-uniform  electric  rotation  can  distort  the  distribution  function 
giving  rise  to  instability  (section  2.1).  The  effects  of  this 
velocity  shear  become  more  pronounced  when  the  electric  rotation 
varies  appreciably  over  a  region  of  the  order  of  the  density  scale 
length.  When  C  is  less  than  zero  (electric  rotation  opposes  the 
diamagnetic  rotation) ,  the  electric  rotation  tends  to  cancel  the  non¬ 
uniformity  in  the  diamagnetic  rotation.  The  eigenfunction  becomes 
delocalized  and  the  growth  rate  drops  (figure  1$.,  equations  3.15  and 

3.16).  When  C„  and  L  have  values  which  best  cancel  the  non- 
E  E 


uniformities  in  the  diamagnetic  rotation,  the  growth  rate  is  minimized 
JR  =  mLNE 


and  WD  =  mCM1_,  (figure  7). 


When 


(W 


i+b  (le/ln) 


<  0,  the  Weber-Hermite  approximation 


becomes  invalid  since  the  "potential"  function  Q  develops  four  turning 


' 
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(cyv 

points.  The  numerical  results  show  that  as  —  becomes 

<VV 

increasingly  negative,  the  growth  rate  again  increases  while  the 
eigenfunction  localizes  over  the  negative  portions  of  the 
"potential  well"  (figure  17) .  This  behavior  in  the  growth  rate 
can  be  understood  from  the  following  argument.  If  equation  (3.1) 
is  multiplied  by  the  complex  conjugate  of  y  and  integrated  from 
0  to  the  edge  of  the  cylinder  r  ,  one  obtains: 


dy 


dr 


dr  + 


*  2 


filn-dr  =  o 


J  4 


(3.18) 


Im(Q) | y |  dr 


0 


o 


z 

By  ignoring  the  small  imaginary  part  of  a^F^^  equation  (3.18)  yields 
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(V^ESW  ly 
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By  letting  only  the  essential  quantities  vary  in  the  region 
where  |y(r) |  is  appreciable,  equation  (3.19)  becomes: 

WI  =  ^  (mHEN  ~  V/  Wp 

\k"v?h 

where  "  "  denotes  an  averaging  process  with  the  eigen¬ 

function  acting  as  a  probability  distribution.  Thus  the  growth 
rate  is  determined  by  what  the  wave  "feels"  over  the  whole  region 
of  the  plasma.  When  Q  has  two  turning  points,  (mW  -W  )  is  always 

Li  IN  O 

positive  over  the  region  of  wave  localization  (figure  13) .  When 
Q  has  four  turning  points,  )  can  become  negative  in  regions 

where  the  wave  amplitude  is  appreciable  (region  1  of  figure  17) 
thus  contributing  a  stabilizing  influence.  However,  the  eigen¬ 
function  is  then  found  to  be  strongly  localized  over  the  most  un¬ 
stable  region  (region  2  of  figure  17) .  This  is  why  a  "potential 
well"  with  four  turning  points  will  not  minimize  the  growth  rate. 

The  growth  rate  is  found  to  be  minimized  when  the  "potential  well" 
is  as  broad  and  shallow  as  possible.  The  analytical  results  (equations 
3.15  and  3.16)  suggest  that  this  occurs  when  the  electric  rotation 
effectively  cancels  the  non-uniformity  in  the  diamagnetic  rotation. 

The  effect  of  a  decrease  in  density  on  the  modes  considered  is 
to  decrease  the  wave  frequency  and  increase  the  growth  rate.  A 
decrease  in  density  seems  to  have  the  most  influence  when  the  electric 
rotation  opposes  the  diamagnetic  rotation  in  such  a  manner  as  to  give 


(3.20) 


' 
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a  minimal  growth  rate.  This  is  best  observed  in  figure  7 

and  figure  9.  However,  decreasing  the  density  even  by  a  factor 

of  ten  does  not  change  the  manner  in  which  the  non-uniform 

rotations  affect  the  growth  rate.  It  would  seem  that  terms 

2 

involving  ( C  / C)  can  be  ignored  with  little  consequence  to  the 

n 

7  8 

qualitative  results  for  densities  in  the  range  5x10  to  5x10 


particles/cm  . 
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CHAPTER  IV 

CONCLUSION 

The  collisionless  drift  wave,  with  a  component  of  the 
propagation  vector  parallel  to  the  magnetic  field,  has  been 
examined  in  a  cylindrical  plasma  for  the  limit  (a^/L^)  <<  1. 

The  differential  equation  describing  the  wave  has  included  the 
effects  of  non-uniformities  in  the  electric  and  diamagnetic 
rotations  as  well  as  a  finite  Alfveh  velocity.  Accurate 
numerical  results  have  been  presented.  Approximate  analytic 
solutions  have  also  been  obtained  using  the  local  approxi¬ 
mation  and  a  Weber-Hermite  approximation.  The  unstable  modes 
found  are  ultimately  due  to  the  finite  Larmor  radius  of  the  ions 
since  when  a.  ->  0  the  wave  becomes  neutrally  stable.  Finite 

l 

Larmor  radius  effects  force  ions  to  experience  the  electric  fields 
and  density  gradients  within  the  plasma  differently  than  the 
electrons  since  the  fields  and  gradients  acting  on  an  ion  must  be 
averaged  over  its  Larmor  radius.  This  causes  the .wave  frequency 
to  be  less  than  the  diamagnetic  frequency  over  all  or  portions  of 
the  wave,  thus  upsetting  the  fine  balance  between  electron  Landau 
damping  and  destabilizing  resonant  particle  effects.  Any  other 
effects  within  the  plasma  which  enhance  the  resonant  particle  effects 
will  destabilize  the  mode,  It  has  been  found  that  non-uniformities 
in  the  electric  and  diamagnetic  rotations  can  provide  the  dominant 
contribution  to  the  growth  rate.  The  effect  of  non-uniform  rotations 


•  • 
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is  further  enhanced  in  lower  density  plasmas  where  terms 

2 

involving  (C  /C)  are  no  longer  negligible.  Thus  instabilities 
occuring  in  the  presence  of  density  gradients  may  in  fact 
derive  most  of  their  energy  from  the  velocity  shear  resulting 
from  non-uniform  rotations.  This  is  thought  to  be  the  case  for 
the  collisional  drift  mode^ .  Since  it  was  shown  that  an 
electric  rotation  opposing  the  diamagnetic  rotation  can  almost 
eliminate  velocity  shear,  imposing  external  radial  electric 
fields  on  a  plasma  column  may  be  an  effective  method  for  stabil¬ 
izing  the  drift  instability. 
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APPENDIX  A 

DERIVATION  OF  THE  DIFFERENTIAL  EQUATION 
DESCRIBING  THE  DRIFT  INSTABILITY 


The  differential  equation  describing  the  drift  instability 
in  a  non-uniform  radial  electric  field  is  derived  below  by 

4 

using  the  guiding  centre  approximation  in  cylindrical  geometry  . 
Ion  Landau  damping,  electron  Landau  damping  and  finite  Larmor 
radius  effects  are  included  in  this  analysis.  All  assumptions 
are  based  on  the  experimental  model  discussed  in  sections  2.2 
and  2.3.  These  assumptions  are  listed  below: 


i.  A  spacially  uniform  axial  magnetic  field  B  =  Be^ 
exists  in  the  plasma. 

ii.  The  velocities  of  the  ions  and  electrons  parallel 
to  the  magnetic  field  are  described  by  Maxwellian 
distribution  functions  at  the  same  uniform 
temperature  T. 

iii .  The  zero  order  radial  electric  field  and  density 
profiles  are  those  given  in  section  2.3. 

iv.  The  electron  diamagnetic  drift  frequency  (W  ) 

and  the  electric  drift  frequency  (Wg)  are  of  the 
same  order  of  magnitude. 

O 

v.  All  terms  of  order  (a^r)  and  higher  are  ignored 


since  the  region  near  the  cylindrical  axis  is  of  no 


. 
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interest  (a  is  the  ion  Larmor  radius  and  r  is  the 
distance  from  the  cylindrical  axis) . 


vi .  A  small  electric  field  perturbation  E^  is  assumed  to 
exist  within  the  plasma  at  time  t  =  0.  E  has  the 
following  properties : 


E^  ^(r)exp  i(m0  +  K,,z  -  Wt) 


(A.  1) 


V  X  E  =  0 


Ei  - 


■V<j). 


i  -»■ 

E, 


<< 


-> 

E 


o 


where 

— y 

E  =  zero  order  radial  electric  field 

o 

m  =  azimuthal  wave  number 

W  =  complex  wave  frequency 

K,,  =  wave  number  parallel  to  the  magnetic  field 

<j)^  =  perturbed  electric  potential 


vii.  The  wave  frequency  divided  by  the  ion  cylotron 

2 

frequency  (Eh)  is  assumed  to  be  of  the  order  (a^/L) 

where  L  is  a  characteristic  length  in  the  plasma 

(density  or  electric  scale  length) .  This  assumption  is 

based  on  the  argument  of  section  2.1  which  shows  that 

3 

W  v  O(W^)  .  Terms  of  order  (a^/L)  and  higher  are 
ignored . 

The  differential  equation  describing  the  instability  is  derived 
in  a  cylindrical  geometry  where  6  ,  £  ,  are  unit  vectors  in  the 


' 

• 

' 
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r,  6  and  z  directions  respectively.  Electromagnetic  units  are 
used  throughout. 

The  zero  order  electric  drift  velocity  perpendicular  to 
the  magnetic  field  is  given  for  electrons  or  ions  as : 


- ^ 
V 


E  X  B 
o 

B2 


dE 


o 


.  B  dt 


1  >e 


(A. 2) 


where  "i,e"  denotes  ions  or  electrons  respectively. 


The  effect  of  a  finite  ion  Larmor  radius  (FLR)  is  obtained  by 


-> 


averaging  the  electric  field  Eq  seen  by  an  ion  over  one  gyration 


period.  This  is  accomplished  (to  the  order  required)  by  applying 
the  operator 


?  2 

1  +  1/ 4a±V 


->  16 
to  E  in  cylindrical  co-ordinates 
o 


Both  the  inertial  and  FLR  effects  of  the  electrons  are  ignored 
since  the  electron  mass  is  much  less  than  the  ionic  mass.  The 
zero  order  drifts  for  electrons  and  ions  are  then  given  by: 


o 


1  +  ajr/4  (V2 


-  l/r2)j  CrWE)fi0 


„2  - 
^E  e0 
ft. 

l 


■  (rWE)€e 


(A.3) 


By  the  use  of  Poisson's  equation,  the  zero  order  electron  density 

(nG'>  is  obtained  in  terms  of  the  zero  order  ion  density  (n1)  . 
o'  ° 


e  l 

n  =  n 
o  o 


1  +  (CA/C)2  1  (r2WE)' 

IT. 

i 


r 


(A. 4) 
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where 


B 


[ 


4ttM.ii 

i 


o] 


1, 

"2 


Alfvdn  velocity 


C  =  velocity  of  light 

=  denotes  differentiation  with  respect  to  r 
The  perturbed  electric  field  lii  gives  rise  to  a  first  order 
drift  velocity  perpendicular  to  B  which  is  given  for  electrons 
or  ions  as : 


-> 

V, 


E  X  B 


4- 


n.  B  dt 


(2o  +  V 


[1] 


(A. 5) 


where  Lb  means  that  only  first  order  terms  in  quantities  such 
as  (a^/r)  ,  (W/fb  )  and  |E^|/|B|  are  retained.  The  reason  that 

E  ,  a  zero  order  term  is  included  in  the  expression  for  V1  is 


o 


that  "d/dt"  is  a  convective  derivative 


3/St  +  (V 


o+vv]- 


Thus  a  term  involving  the  zero  order  electric  field  when 
multiplied  by  a  first  order  velocity  will  give  rise  to  a  term 
which  is  of  the  first  order.  Equation  A. 5  is  evaluated  using 
the  finite  Larmor  radius  operator  for  the  ions.  In  calculating 
the  first  order  electron  drift  velocity,  FLR  and  inertial  effects 


are  ignored.  The  following  expressions  are  then  obtained  for 


the  first  order  ion  and  electron  drifts: 


■ 


■ 
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1  -  2W£  +  a2(V2-l/r2)  1 

~  4~" 

1 


~y 

E  X  B 


B 


O  O 

i  (W  -  mWI  -  mafsi./2r  )E, 
STB  l  a  1 

1 


£r  r(WE/fll),Eie/B 


B 


X  B 
2 


where 


WI 


V^VV 


+  a2/4(W„"  +  3W  ') 
X  x  —  X 
r 


v2r>  =  e  72(  >  +  v2(  >8  +  ezv2(  >z 


(A.  6) 


In  order  to  calculate  the  first  order  particle  density,  an 
expression  for  the  first  order  guiding  centre  density  is 
required.  This  is  obtained  by  using  the  equation  of  continuity 
for  guiding  centres. 


a_ 

at 


(N  +N )  +  (V  +V1 ) • V (N  +N,)  +  (N  +N, ) V* (V  +V, )  =  0 
v  o  1  o  1  o  1  o  1  o  1 


(A.  7) 


Using  the  assumptions  that: 
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x . 


varies  as  exp  i(m0-K,,z-Wt) 


•  • 

11.  N  and  V  are  functions  of  r  only  with  V  in 
o  o  J  o 


the  i. 


0 


in . 


Second,  order  terms  V^'VN^,  can  be  ignored 


the  equation  of  continuity  becomes: 


•iWN1  +  (V  )nlmNl  +  (V,)  N*  +  N  V.V,  =  0 

1  o  0 - -  1  r  o  o  1 


(A.  8) 


When  equations  (A. 3)  and  A. 6)  are  substituted  into  equation 
(A. 8),  the  following  expressions  for  the  perturbed  ion  and 
electron  guiding  centre  densities  are  obtained: 


N. 


-mN' 


o 


(W-mWI)rB 


2  2 

2W„  a .  Vx 

— I  +  a _ 

Si,  4 


<J> ,  Vx  *  N  V  <£ 

1  ,  o  1 

r  fi.B 

l 


i  '  i  »  ' 

3mN  (J)1  W„  imcj)  (N  W  ) 

O  1  B  .  1  O  B 

+  -7TZ — - — rrx - r 


iN1  d(vh 
o  1  z 


(W-mW  )  rB .  ( W-mWI )  S! .  B 

B  1  l 


(W-mWI)  dz 


N, 


-mN 


iN 


o 


o 


d(vIh 


(W-mW  ) rB 

ht 


(W-mWg)  dz 


(A. 9) 


(Terras  of  the  order  (M  /M  )  have  been  ignored.) 


From  Schmid tXD  the  particle  density  can  be  expressed  in  terms  of 
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the  guiding  centre  density  by  applying  the  operator 


2  ? 

1  +  l/4af V 


to  .  The  perturbed  ion  and  electron  particle  densities 
obtained  from  equation  (A. 9)  are  then  given  by: 


en„ 


2  2  2  • 

e  J  1  d  Sd^  (1-m  )Sip  rW  n^ 

M  n  2  ' r  dr  dr  3 

M .  .  W  r  r 

l  l  o 


mW  fi.rn1^  -  ifi.Brn1  d(V^)  , 
O  1  O  l  o  1  z 


dz 


e 

en . 
i 


-e 


M.fl.  W  r 
11  o 


y  rW  ft.mn6  +  ifi.Brn6  d(V*f) 

<  o  1  °  1  o  1  z 


dz 


J 


(A. 10) 


where 

e  =  electron  charge 

W  =  W  -  mW„ 

o  E 

V  =  <Pi/W  r 

1  o 

O  O  * 

S  =  W  r  n(l  “  mWjW  ) 
o  o  o 

=  ion  diamagnetic  drift  frequency  (equation  2.24) 
W  =  electric  drift  frequency  (equation  2.23) 

ill 


The  final  differential  equation  will  be  expressed  in  terms  of  the 
variable  \p  for  which  physical  interpretation  will  now  be  given. 
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The  equation  of  motion  for  the  guiding  centres,  ignoring 
inertial  drifts,  is  given  by: 


dr 
_ £ 

dt 


yE 

B 


r  -  r  +  r  and  E  -  E  +  E1  where  r  ..  is  the  displacement 
g  g°  gl  o  1  gl  r 

of  a  guiding  centre  in  a  reference  frame  moving  at  the  zero  order 

drift  velocity.  Then  the  first  order  equation  is: 


(3/3t  +  V  • V) r  . 

o  gl 


vq 

—  (Ei  "  -v*i> 


Using  the  approximation  that  r  ,  «  enexp  i(m0-Wt) 

gl  0 


the  0-component 


of  this  equation  gives: 


-m 

B 


W  r 


o 


W  r 


o 


Thus  is  proportional  to  the  perturbed  guiding  centre  displacement 
resulting  from  the  perturbation  <f>^. 

The  terms  which  still  remain  to  be  evaluated  in  equation  (A. 10) 
d<^)z  d(V^)z 

are  — -  and  — - .  In  this  analysis  it  is  assumed  that  the 

dz  dz 

electron-ion  collision  frequency  v  <<  fh  ,  and  that  the  ion 
Larmor  radius  is  much  smaller  than  the  wavelength  of  the  perturbation. 
Therefore,  the  drift  approximation  to  the  Vlasov  equation 


can  be 


67 


d(V  ) 

-L  z 

used  to  calculate  — - .  This  essentially  replaces  the  particles 

by  quasi-particles  located  at  the  guiding  centres  with  magnetic 

moments  y.  The  distribution  function  f  is  now  considered  to  be  a 

•> 

function  of  the  guiding  centre  position  r  ,  the  parallel  velocity 

§ 

and  the  magnetic  moment  y.  It  is  not  considered  a  function  of 

the  guiding  centre  velocity  since  the  guiding  centre  velocity  is 

determined  by  r  (equation  A. 12).  In  "drift-space",  the  Vlasov 

g 

equation  becomes: 


M  =  M  +  ji  Ik  +  it. 

dt  9t  g  dt  9 V  dt  9y  dt 

z 

where  f  =  f(r  ,V  ,y) . 

K  2 


o 


(A. 11) 


The  equations  of  motion  for  the  guiding  centres  are  given  by: 


dr  (1  xl)  xS 

g  =  _ z  _o _ 1 _ 

dt  "  B  g2 


dV 

z 

dt 


^i>» 

M. 

i,e 


e  for  ions 
-e  for  electrons 


(A. 12) 


—  =  0  (consistent  with  the  guiding  centre  approximation) 
dt 


Equation  (A. 11)  is  linearized  by  assuming  f  =  fQ  +  f±  where  ^  is 


■ 
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a  small  first  order  perturbation  and  f  is  a  zero  order 

o 

Maxwellian  given  by: 


-  —r 

M. 

H 

-M.  V  2' 

i,e 

exp 

l ,  e  z 

27TkT 

2kT 

L  J 

By  substituting  equations  (A. 12)  into  (A. 11)  the  following 
expression  is  obtained  for  f^: 


lq(Vz 

kT 


mW. 


V  - 


D 


z  K„ 


W  -  mW  -  K„V 
E  z 


o 


(A. 13) 


where 


n 


W. 


o  -kT 


D 


n  r  qB 
o 


(electron  or  ion  diamagnetic  frequency) 


The  macroscopic  velocity  parallel  to  B  is  obtained  by  averaging 

V  over  velocity  space,  that  is: 
z 


(h>z 


oo 


1_ 

n 

o 


f.V  dV 
1  z  z 


CO 


(Methods  to  evaluate  this  integral 
are  given  in  reference  17) 


>1, 


Thus  the  expressions  obtained  for  n  d(V..)  and  n  d(V,)  are: 
r  o  1  z  o  1  z 


dz 


dz 
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ne  d(V^) 
o  r  z 


dz 


-2efVA 

MeV™2  K" 


1  +  i/ir£  E(£  ) 

s:  » 

1 

g 

1 

e  e 

O  EH 

no  d(^i) 

°  1  z 

dz 


+2e(E  )  n 
1  z  o 

mi42  K" 


1  +  i/ir£.E(£.) 

i  i 


W  -  mW , 
o 


(A. 14) 


where 


VTU  =  2kT/M 


TH 


v1 

TH 


k 


W. 


EN 


E(0  = 


2kT/M± 


Wo/K"VTH 


Vk"4 


(electron  thermal  velocity) 
(.ion  thermal  velocity) 


electron  diamagnetic  frequency  (Equation  2.26) 

15 


2 

Ttt 


2  r2 

~S  A 

e  ds  e 


«  OO 


The  final  equation  necessary  for  a  self-consistent  set  is 

Poisson's  equation. 

_2  ,  .  ^,2,  i  e. 

-V  cfi  =  4ttC  (n^  -  n^)e 


(A. 15) 
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By  substituting  the  results  of  equations  (A.l),  (A. 10)  and 
(A. 14)  into  (A.  .15)  ,  the  following  equation  is  obtained  for  \p: 


!  \ 

Id  IS  df  +  \p 
7  dr  i^°  drj 


(l-m2)S  +  rW2nX 

: - o -  O  o 


2  2  i  ?  2 
(C./C)  K„nVr 
A  o  o 


-  ip 


( 


2r\  K(Vm,V  z(Wo/K"vra)  +  no<VmW~>  z(wo/K"vra> 


2  \ 

a. 

1 


i  A 

H 

/ 


(A.  16) 


where 


o 


2  3  i 
W  r  n 
o  o 


1  +  (CA/C)  -  mW„/Wo 


z(0 


l  +  i  Ace  (5)1 
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APPENDIX  B 

NUMERICAL  METHOD  OF  SOLUTION 

In  this  section  the  numerical  method  of  solution  for 
the  equation  discussed  in  chapter  three  will  be  given.  The 
equation  is  of  the  form: 

y"(r)  -  Q(r,W)y(r)  =  0-  (B.l) 

where  Q  is  an  arbitrary  function  of  r  and  the  complex  eigen¬ 
value  W.  The  assumptions  that  the  mode  is  localized  and  that 
a  perfectly  conducting  boundary  exists  at  r  =  r  (plasma  radius) 
gives  rise  to  the  following  boundary  conditions: 
i.  y (0)  =  0 

ii.  y(r  )  -  0 

P 

The  finite  difference  form  of  equation  (B.l)  is  given  by: 

yn+i  _  (2  +  h\(w))  +  yn-i  =  0  (b.2) 

where 

h  =  step  size 

yn  =  with  n  an  integer 

Qn  =  Q(nh) 

If  the  variable  N  is  taken  to  represent  the  number  of  internal 
node  points,  then  N  equations  of  the  form  given  by  (B.2)  may  be 


written : 
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y* 


(2+h  Q1)y1 


(2+h  Q2)y2 


+ 


y 


o 


0 

0 


N+l 


(2+h  VyN 


+  y 


N-l 


0 


with  yo  =  yN+1  =  0 


The  above  equations  may  be  expressed  in  matrix  form  as: 


M(W)  V  =  0 


(B .  3) 


wnere 


M(W)  = 


(2+h  Q1(W))  1 


0 


0 


-(2+h  Q2(W)) 


0 


\ 


0 

f 


(2+h  Qn_1(W)) 


_  1 


0 

» 


0 


-(2+h  Qn(W)) 


V  = 
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A  method  for  solving  this  non-linear  eigenvalue  problem  has 


been  derived  by  M.  Osborne 


13 


The  method  has  the  following 


assumptions  and  properties: 

1.  The  method  is  second  order  in  that  the  error  <5 .  , , 

1+1 

2 

at  the  j+1  stage  is  proportional  to  6.  where  6. 

1  1 

,  .  .  th 

is  the  error  at  the  j  stage. 

2.  The  eigenvalues  are  assumed  to  be  distinct. 

3.  Initial  approximations  are  required.  If 

V.  .  .  ,  =  V  +  q  where  llqll  <6 

—initial  —true  -1  1 


W.  .  .  ,  =  W  +  <$n 

initial  true 


then  the  conditions  for  convergence  are  that: 
i .  8  <  1/3 

ii.  ] n |  <  1.5a  where  a  is  the  minimum 
separation  of  the  eigenvalues. 

4.  The  iteration  scheme  used  is  given  by: 


M(W.)  V.,,  dM(W.)  V 

l  —i+l  =  i—i 


dW 


(V.) 

-i  p 


(B .  4) 


W. , 1  =  W.  -  1/(V  ) 

i+l  l  —i+l  p 


where  p  is  the  index  of  the  component  of  maximum 


modulus  in  V. 


. 


. 


» 

- 
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A  computer  program  utilizing  the  above  numerical  method 
was  written  in  Fortran  IV.  The  initial  eigenvectors  were 
assumed  to  have  components  of  the  following  form: 


2 ,  N  2r  . 


V  =  r  (r-r  ) 
k  p 


sin(nTrr/r  )  +  i  sin(  (n+1)  irr/r  ) 
P  P 


where 


r  = 
P 


kh 

(N+l)h 


n  =  radial  mode  number 


k  =  component  index  of  V  (k=l,2...N) 


The  initial  values  for  W  were  chosen  as  Max(mI7  T  +  mlJ )  .  These 

EN  E 

initial  values  were  systematically  reduced  until  one  of  them 
converged  to  the  desired  radial  mode.  The  frequencies,  growth 
rates  and  eigenvectors  were  all  calculated  with  a  relative  error 


less  than  .0002. 


■ 


